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During this period we initiated work on modeled PDF, i.e., PDFs that are governed by an 
evolution equation. An abstract 1 was submitted to the AIAA/SAM/ASME/ASEE 30th Joint 
Propulsion Conference and Exhibit. A copy of the abstract is enclosed. 


1. Baurle, R. A., Hsu, A. T., and Hassan, H. A., "Comparison of Assumed and Evolution PDF' 
in Supersonic Turbulent Combustion Calculations", Abstract submitted to the 30th Joint 
Propulsion Conference and Exhibit. 
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The objective of this study is to compare the use of assumed pdf (probability density 
function) approaches for modeling supersonic turbulent reacting flowfields with the more 
elaborate approach where the pdf evolution equation is solved. Assumed pdf approaches 
for averaging the chemical source terms require modest increases in CPU time typically of 
the order of 20 % above treating the source terms as “laminar”. However, it is difficult 
to assume a form for these pdf’s a priori that correctly mimics the behavior of the actual 
pdf governing the flow. Solving the evolution equation for the pdf is a theoretically sound 
approach, but because of the large dimensionality of this function, its solution requires a 
Monte Carlo method which is computationally expensive and slow to converge. Preliminary 
results show both pdf approaches to yield similar solutions for the mean flow variables. 


Introduction 

In order to obtain accurate predictions for turbulent reacting flowfields, the effect of 
turbulence on the chemical source terms must be addressed. These source terms are highly 
nonlinear functions of both temperature and composition and as a result employing typical 
time averaging is difficult. An attractive alternative for obtaining these averages is to utilize 
pdf models which can be used to obtain averages for functions of any form. These pdf 
approaches consist of two types. The first being the assumed pdf approach where the pdf 
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is completely defined knowing only the first couple of moments, and the second approach 
where an evolution equation is solved for the pdf. 

Many assumed forms for pdf’s have been used for obtaining averages for the chemical 
source terms in supersonic reacting flows. All of these forms have either neglected the effects 
of composition fluctuations, or assumed statistical independence between the temperature 
and the composition. Frankel et al. 1 examined the effects of temperature fluctuations on 
the chemical source terms using both Gaussian and Beta forms for the pdf. Both assumed 
forms yielded similar results and showed modest improvements over treating these source 
terms in a laminar fashion when compared to experiment. Later Narayan and Girimaji 2 
combined a moment method to account for temperature fluctuations with a multi-variate 
Beta distribution developed by Girimaji 3 to account for composition fluctuations. Their 
results showed the effect of the multi-variate Beta distribution to be minimal on the mean 
flow variables. This same behavior was seen later by Baurle et al. 4 

The use of evolution pdf methods for low speed turbulent reacting flows have been exam- 
ined extensively. 5 ’ 6 ’ 7 However, until recently 8 , the extension of these methods to high speed 
flows has been practically non-existent due to difficulties arising from shocks and strong di- 
latation terms. Hsu et al. 8 solved the evolution equation for the joint pdf of specific enthalpy 
and mass fractions in conjunction with a CFD flow solver for the Navier-Stokes equations. 
Significant improvements over calculations without the pdf were noted when compared to 
experimental data for high speed reacting flows. 


The Assumed Joint PDF Approach 


For the assumed pdf approach, the species continuity equations are solved in conjunction 
with the Navier-Stokes equations. To obtain an average of the chemical source terms (d>k) 
that appear in the species equations, 


Uk = 






( 1 ) 


a joint pdf of temperature and composition is needed. In this work, the joint pdf was chosen 
as 

V(T,Y k )=V 1 (T)V 2 (Y k ) (2) 

Here, the pdf of temperature ('Pi(T)) was chosen as a Gaussian distribution 

(T-f) 2 


V 1 (T) = 


\J 2tT"‘‘ 


■ exp 


2T" 2 


( 3 ) 


and the pdf of the mass fractions (^(V*)) was chosen as the multi-variate Beta distribution 
developed by Girimaji 3 
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It can be shown that 
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To use this approach the solution of two extra transport equations are required. One for 
the enthalpy (or temperature) variance and one for the species mass fraction variance sum 

M- 


The Evolution Joint PDF Approach 

For the evolution pdf approach, the species continuity equations are replaced by the pdf 
evolution equation. The simulation of this equation yields the mean mass fractions which is 
all that is required by the Reynold’s averaged Navier-Stokes equations. The single point joint 
pdf evolution equation is derived from the specific enthalpy and species continuity equations. 
If the dissipation due to viscosity is neglected, this equation can be written as 


d[pV ) d(pujV ) d(io k V} 
dt ^ dxj dipk 


d{p(u';\ i, k ,6)v) a((Sf \MP) 
dxi do 


(7) 


where M (il>k,0\ x,t) represents the effect of molecular mixing, and (if>k,0) are the sample 
space variables corresponding to the random variables ( Y k ,h ). The terms on the right hand 
side of the above equation require modeling which is discussed in references 7, 8, and 9. 

Due to the large dimensionality of the joint pdf, solving this equation with a finite differ- 
ence approach is not practical since the computational cost can be shown to rise exponentially 
with the dimensionality of the pdf. The computational cost of a Monte Carlo simulation of 
the above expression can be shown to rise linearly with the dimensionality of the pdf. 10 Thus, 
the above equation is simulated using a fractional step Monte Carlo scheme as described in 
reference 10. 


Results and Discussion 

The cases considered involved the parallel injection of hydrogen into a supersonic stream 
of air as shown by the schematic in Figure 1, and results are compared with the experimental 
data of Beach et al. 11 Table 1 gives a summary of the burner exit conditions. Measurements 
were taken at axial stations of 8.26, 15.5, 21.7, and 27.9 inner diameters downstream of the 
burner. All computations were carried out using a 61 X 71 grid on a domain of 30 X 2 inner 
diameters. 

The CFD solver employs a cell centered finite volume approach stepped in time using a 
4 stage Runge-Kutta scheme. To alleviate the stiffness caused by the different time scales 
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between the fluid motion and the chemistry, the chemical source terms are treated in a point 
implicit manner. 

For the assumed pdf approach, the 7 step kinetic model by Jachimowski 12 was employed 
as given in Table 2. Figure 2 compares the calculated mean mass fraction profiles with the 
experimental data at all 4 axial stations where measurements were taken. As can be seen, 
fair agreement with experiment is noted at each station. The total cpu time to reach the 
steady state (5000 iterations) was approximately 860 seconds on a Cray C90. 

For the evolution pdf approach, a simpler 2 step kinetic model by Rogers and Chinitz 13 
was employed as given in Table 3. Once again good agreement with experiment is noted, in 
fact the results are very close to that predicted by the assumed pdf approach where a more 
elaborate chemistry model was used. One noticeable difference is the location of the peak 
in the H 2 0 profile comparisons. It is shifted down somewhat more than in the assumed pdf 
solution which is consistent with the experimental measurements. The total cpu time to 
reach the steady state (6000 iterations) was approximately 2640 seconds on a Cray C90. 

The final paper will implement a general chemistry model into the pdf Monte Carlo solver 
so a more meaningful comparison between computational costs and accuracy can be made. 
Also, more detailed comparisons with the above experiment and comparisons with a more 
recent experiment conducted by Cheng et al. 14 will be presented. 
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Tables 


Table 1: Burner Exit Conditions 


Exit Conditions 

Inner Jet 

Outer Jet 

Mach Number 

2.0 

1.9 

Temperature [°K] 

251.0 

1495.0 

Pressure [MPa] 

0.1 

0.1 

Mass Fraction 



Yh 2 

1.0 

0.0 

Yo 2 

0.0 

0.241 

Yn 2 

0.0 

0.478 

Yh 2 o 

0.0 

0.281 

all others 

0.0 

0.0 


Inner Jet Diameter = 0.006525 m 
Lip Thickness = 0.0015 m 
Outer Jet Diameter = 0.0653 m 


Table 2: Abridged Jachimowski Chemistry Model 


Reaction 

A 

b 

T a 

i?2 "H O 2 ^ OH T OH 

1.70E+13 

0.0 

24157.0 

H + 0 2 - OH + 0 

1.20E+17 

-0.91 

8310.5 

OH T H 2 ^ H 2 O T H 

2.20E+13 

0.0 

2591.8 

O + H 2 ^ OH + H 

5.06E+04 

2.67 

3165.6 

OH + OH ^ H 2 0 + 0 

6.30E+12 

0.0 

548.6 

H + OH+M ¥* H 2 0 + M 

2.21E+22 

-2.0 

0.0 

H + H+M^H 2 + M 

7.30E+17 

-1.0 

0.0 


Units of A are a multiple of m ^. sec 


T able 3: Rogers and Chinitz Chemistry Mod el 


Reaction 

A 

b 

T 

H 2 + 0 2 ^ 2 OH 

Ax 

-10.0 

2448.0 

H 2 + 2 OH ^ 2 H 2 0 

a 2 

-13.0 

21389.0 


Units of A are a multiple of m ^, sec 
Ax = + 8.917^ - 28.95) 10 47 

A 2 = - 0.833<£ + 2.0) 10 64 
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Figures 



Figure 1 . Schematic of Test Apparatus 


7 



' 0.00 0.10 0.20 0.30 0.40 0 50 0.60 0.70 0.80 0.90 1.00 ' 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0 90 1 00 



